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The goal of this Section is to formulate some of the basic results on 
the theory of integral equations and mention some of its applications. The 
literature of this subject is very large. Proofs are not given due to the space 
restriction. The results are taken from the works mentioned in the references. 


1. Fredholm equations 


1.1. Fredholm alternative 

One of the most important results of the theory of integral equations is 
the Fredholm alternative. The results in this Subsection are taken from 151, 


20 , 21 


Consider the equation 


u = Ku + /, Ku = / Ii(x,y)u(y)dy. (1) 

J D 

Here / and K ( x , y) are given functions, K(x, y) is called the kernel of the 
operator Ii. Equation (□]) is considered usually in H = L 2 (D) or C(D). 

1 /Q 

The first is a Hilbert space with the norm ||u|| := ff D \u\ 2 dx) and the 
second is a Banach space with the norm ||w|| = sup xeD |u(a;)|. If the domain 
D C W' is assumed bounded, then K is compact in H if, for example, 
f D f D | K(x, y)\ 2 dxdy < oo, and K is compact in C(D) if, for example, Ii(x, y) 
is a continuous function, or, more generally, sup x£ D f D \K(x,y)\dy < oo. 

For equation (Jl]) with compact operators, the basic result is the Fredholm 
alternative. To formulate it one needs some notations. Consider equation 
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© with a parameter p e C: 


u = fiKu + /, (2) 

and its adjoint equation 

v = JlK*v + g, K*(x,y ) = K(y,x), (3) 

where the overbar stands for complex conjugate. By (2 0 ) and (3o) the cor¬ 
responding equations with / = 0, respectively, g — 0 are denoted. If K is 
compact so is K*. In H one has ( Ku,v ) = ( u,K*v ), where (u,v) = j D uvdx. 
If equation (2 0 ) has a non-trivial solution u 0 , u 0 ^ 0, then the corresponding 
parameter p is called the characteristic value of K and the non-trivial solu¬ 
tion «o is called an eigenfunction of K. The value A = ^ is then called an 
eigenvalue of K. Similar terminology is used for equation (3 0 ). Let us denote 
by N(I — fiK) = {u : u = gKu }, R(I — p/L) = {/:/ = (/ — gK)u}. Let us 
formulate the Fredholm alternative: 


Theorem 1.1. If N(I - /iK) = {0}, then N(I - ftK*) = {0}, R(I - fiK) = 
H, R(I - jlK*) = H. 

If N (/ — ixK) 7 ^ {0}, then dim N(I — /xK) = n < oo, dim A r (I — JlK*) = 
n, f € R(I — /iK) if and only if f _L N(I — JlK*), and g e R(I — JlK*) if 
and only if g T N(I — fiK). 


Remark 1.2. Conditions / _!_ N(I — JlK*) and g _L N(I — fiK) are often 
written as ( f,v 0 j ) = 0, 1 < j < n, and, respectively, ( g,u 0 j) — 0, 1 < j < n, 
where {;%}” =1 is a basis of the subspace N (/ — JlK *), and {woi}j=i is a basis 
of the subspace N(I — gK). 


Remark 1.3. Theorem 1 1.1 1 is very useful in applications. It says that equation 
dU) is solvable for any / as long as equation (2q) has only the trivial solution 


t>o = 0. A simple proof of Theorem II.II can be found in [20] and 21 


Remark 1.4. Theorem 11.11 can be stated similarly for equations (] 2 ]) and (J3J) 
in C{D). 


Remark 1.5. If n > 0 and / G R(I — gK), then the general solution to 

n 

equation ([ 2 ]) is u — u + CjU 0 j, where Cj, 1 < j < n, are arbitrary constants 

3 = 1 

and u is a particular solution to (J 2 ]). Similar result holds for v. 
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1.2. Degenerate kernels 
Suppose that 


M 

K(x,y) = y^a m (x)b m (y), (4) 

m= 1 

where the functions a m and b m are linearly independent, a m , b m G H. The 
operator K with the degenerate kernel (|D) is called a finite-rank operator. 
Equation ([2]) with degenerate kernel (JJJ) can be reduced to a linear algebraic 
system. Denote 

u m = ub m dx. (5) 

J D 

Multiply equation (J2]) by b p , integrate over D, and get 

M 

Hp R ^ ^ (, 'pm ^ru T fpi 1 ^ P ^ Tf, (6) 

771=1 

where a pm := ( a m ,b p ). There is a one-to-one correspondence between solu¬ 
tions to equation (J2J) and (151) : 

M 

u(x) = n a m (x)u m + f(x). (7) 

771=1 


Exercise 1. Solve the equation u{x) = fi / sin(x — t)u(t)dt + 1. 

Jo 

Hint: Use the formula sin (x — t) = sin x cos t, — cos x sin t. Find charac¬ 
teristic values and eigenfunctions of the operator with the kernel sin(x — t) 
on the interval ( 0 , 7 r). 

Remark 1.6. If K is a compact integral operator in L 2 (D) then there is a 
sequence K n (x,y) of degenerate kernels such that lim \\K — K n \\ = 0, where 

n—yoo 


K || is the norm of the operator K, ||/l|| := sup 

u^O 


\\Ku\[ 

I Ml 


1.3. Volterra equations 

If K(x,y) = 0 for y > x and dim D = 1 then equation (J2J) is called 
Volterra equation. It is of the form 

PX 

u{x) = n / K(x,y)u(y)dy + f(x), a<x<b. ( 8 ) 

J a 


3 




This is a particular case of the Fredholm equations. But the Volterra equation 
USD h as special property. 

Theorem 1.7. Equation 0) is uniquely solvable for any /i and any f in 
C(D), if the function K(x,y ) is continuous in the region a < y < x, a < 
x < b. The solution to 0) can be obtained by iterations: 

pX 

u n+1 (x) = pVu n + f, ui = /; Vu: = K(x,y)u(y)dy. (9) 

J a 

1-4■ Self adjoint operators 

If K* = K then K is called selfadjoint. In this case equation (0) has at 
least one characteristic value, all characteristic values are real numbers, and 
the corresponding eigenfunctions are orthogonal. One can solve equation (J2]) 
with selfadjoint operator K by the formula 


u = 


OO r. 

,=i Tj ~ A 


(/> u 0j)i 


( 10 ) 


where Uoj = HjKu 0 j , j = 1, 2,.... 

For the eigenvalues of a selfadjoint compact operator K the following result 
(minimax representation) holds. We write K > 0 if (Ku,u) > 0,Vm 6 //, 
and A < B if ( Au,u ) < (£>«, w), Vu G H. 

Theorem 1.8. If K = K* > 0 is compact in H, Kuj = A jUj, Uj ^ 0 then 

(. Ku,u) . , , 

Aj+i = max , j = 0,1,... (11) 

u±Lj \\U\\ Z 

Here L 0 = {0}, Lj is the subspace with the basis {wi,..., Uj}, ||wjj| = 1, (iq, u m ) 
film- Also 

, . (Ku,u) . . . 

A;+i = nun max—-———, dim M — i, 12 

J M u±M \\u\\ 2 

where M runs through the set of j-dimensional subspaces of H. 

One has Ai > A 2 > • • • > 0. The set of all eigenfunctions of K , including 
the ones corresponding to the zero eigenvalue, if such an eigenvalue exists, 
forms an orthonormal basis of H. 

Remark 1.9. If 0 < A < B are compact operators in H, then Xj(A) < A j(B). 
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Remark 1.10. If A is a linear compact operator in H , then A = I/4|, where 
\A\ = {A*A) 1 / 2 > 0 and U is a partial isometry, 

U:R(\A\)^ R(A), \\U\A\f\\ = H4/1I V/ G H, N(U) = N(A ), 

and i?(4l) = R(A*), where the overline denotes the closure in H. One has 

r(A) 

Au = E s j( A )( u , u j) v j- ( 13 ) 

3 = 1 

Here r(A) = dim.R(4l)> Sj > 0 are the eigenvalues of 41- They are called 
the s-values of A, {uj} and {vj} are orthonormal systems, \A\uj = SjU 3 , 
Si > S2 > ■ ■ ■ , Vj = UUj. 

Remark 1.11. If A is a linear compact operator in H then 

Sj+i( A ) = min |4-A'||, (14) 

KGKj 

where K 3 is the set of finite rank operators of dimension < j, that is, 
dimR(Kj) = j. 


1.5. Equations of the first kind 

If Ku = Xu + f and A = 0, then one has 


Ku = f. 


(15) 


If K is compact in II. then equation (ITS]) is called Fredholm equation of the 
first kind. The linear set R(K ) is not closed unless K is a finite-rank operator. 
Therefore, small pertubation of / G R{K) may lead to an equation which 
has no solutions, or, if it has a solution, this solution may differ very much 



. The usual statement of the problem of solving an ill-posed equation (1X5)1 
consists of the following. It is assumed that equation (TT5]l is solvable for the 
exact data /, possibly non-uniquely, that the exact data / are not known 
but the “noisy” data f$ are known, and \\fs — f\\ < S, where <5 > 0 is known. 
Given f$, one wants to find u$ such that lim ||w — u^|| = 0. 

5—>0 


There are many methods to solve the above problem: the DSM (Dynam¬ 
ical Systems Method) [l7| 18], iterative methods, variational regularization, 


and their variations, see [4] and 15 
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Examples of ill-posed problems of practical interest include stable nu¬ 
merical differentiation, stable summation of the Fourier series and integrals 
with perturbed coefficients, stable solution to linear algebraic systems with 
large condition numbers, solving Fredholm and Volterra integral equations 
of the first kind, various deconvolution problems, various inverse problems, 
the Cauchy problem for the Laplace equation and other elliptic equations, 
the backward heat problem, tomography, and many other problems, see flil ]. 

1.6. Numerical solution of Fredholm integral equations 

Let us describe the projection method for solving equation (J2]) with com¬ 
pact operator K. Let L n C H be a linear subspace of H, dim L n = n. 
Assume throughout that the sequence L n C L n +i is limit dense in H , that 
is, lirn p(/, L n ) = 0 for any / G H, where p(f,L n ) = inf \\f - u\\. Let P n 

a—>-00 u(fL n 

be an orthogonal projection operator on L n , that is, P% = P nj Pnf £ L n , 
11/ - Pnf || < ||/ - u\\,Vu in L n . 

Assume that N{I — K ) = {0}. Then ||(J — Jt) _1 || < c. Consider the 
projection method for solving equation ([Tj) : 

P n (I - K)P n u n = P n f. (16) 

Theorem 1.12. If N(I — K) = {0}, K is compact, and {L n } is limit 
dense, then equation m is uniquely solvable for all n sufficiently large, 
and lim \\u n — w|| = 0, where u is the solution to equation W- 

n —^00 

Remark 1.13. Projection method is valid also for solving equations Au = /, 
where A : X —> Y is a linear operator, X and Y are Hilbert spaces. Let 
Y n c Y and X n C A" be n-dimensional subspaces, P n and Q n are projection 
operators on X n and Y n respectively. Assume that the sequences { X n } and 
{Y n } are limit dense and consider the projection equation 

QnAP n U n — Qnfi 'U-n £ A n . (17) 

Theorem 1.14. Assume that A is a bounded linear operator, ||A _1 || < c, 
and 

\\Q n AP n u\\ > c\\P n u\\, VuG A, (18) 

where c > 0 stands for various positive constants independent of n and u. 
Then equation 0 is uniquely solvable for any f e Y for all sufficiently 
large n, and 

lim ||u — u n \| = 0. (19) 

n—>oo 
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Conversely, if equation is uniquely solvable for any f G Y for all suffi¬ 
ciently large n, and HOI) holds, then ITSl) holds. 

Remark 1.15. If (fTBD holds for A, then it will hold for A + B if ||i?|| > 0 is 
sufficiently small. If (TT5jl holds for A, then it holds for A + B if B is compact 
and 11(A + T?) _1 || < c. 

Iterative methods for solving integral equation (jT|) are popular. If ||A'|| < 
1, then equation (JT() has a unique solution for every / G H, and this solution 
can be obtained by the iterative method 


u n+ 1 = Ku n + /, lim 11 u n — ii|| =0, 

n—> oo 


( 20 ) 


where U\ G H is arbitrary. 

Define the spectral radius of a linear bounded operator A in a Banach 
space by the formula p(A) = lim ||24 n || 1 ^ n . It is known that this limit exists. 

n—> oo 

If a is a number, then p(aA) = \a\p(A). The basic fact is: 

If |A| > p(A), then the equation Xu = Au + f is uniquely solvable for any 
f, and the iterative process Xu n+ i = Au n + f converges to u = (A/ — A) -1 / 
for any f G H and any initial approximation u\. 

Given an equation Au = /, can one transform it into an equation u = 
Bu + g with p(B) < 1, so that it can be solved by iterations? To do so, 
one may look for an operator T such that ( A + T)u = Tu + /, (A + T)~ l is 
bounded, so u — (A + T) _1 Tn+ (A + T)" 1 /, and p((A + T) _1 T) < 1. Finding 
such an operator T is, in general, not simple. If A = A*, 0 < m < A < M, 
where m and M are constants, one may transform equation Au = f to the 


2 Au 


+ 


r 7, and get p(I 


2 A 


< 


M—m 


< 1 . 


equivalent form u = u m+M , m+M , , — m+M > - M+m 

If A A* then equation A*Au = A* f is useful since A*A is a sclfadjoint 
operator. The above results can be found, for example, in [§| and [ 27 J. 

Equation Au = f is often convenient to consider in a cone of a Banach 


space A", see 27]. A cone is a closed convex set K C A" which contains with 
a point u all the points tu,t > 0, and such that u G K and — u G K imply 
u = 0. One writes u < v if v — u G K. This defines a semiorder in X. If K 
is invariant with respect to K, that is, u G K implies Au G K, then one can 
give conditions for the solvability of the equation Au = f in K. For example, 
let X = C(D), D C R m is a bounded domain, K = {u : u > 0}, 


Au 


A(x,y)u(y)dy, 


'D 


A(x,y ) > 0. 
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The following result is useful: 

Suppose that there exist functions v and w in K such that Av > v, Aw < 
w, and v < w. Let us assume that every monotone bounded sequence in K, 
ui < U 2 < ■ ■ ■ < U has a limit: 3 u 6 K such that lim \\u n — u\\ = 0. Then 


there exists u G K such that u = Au, where v < u < w. 

A cone is called normal if 0 < u < v implies ||u|| < c||v||, where c > 0 is 
constant independent of u and v. 

2. Integral equations with special kernels 

2.1. Equations with displacement kernel 

The results in this Subsection are taken from [3], [ 27 J. 

Let us mention some equations with displacement kernels. We start with 
the equation 


A u(t) — K(t — s)u(s)ds — f(t), t> 0. 



This equation is solved by taking Laplace transform: Xu — Ku = /, u 
(A — K ) -1 /, where u := Lu := / 0 °° e~ pt u(t)dt, and 



pcT+ioo 
J (T —ZOO 



Equation 



OO 


K(t — s)u(s)ds = /(f), — 00 < t < 00 , 


— OO 


is solved by taking the Fourier transform: u(f) := jf° u(t)e l &dt. One has 
*(0 = (A - K(£))~ l fi and 



In these formal calculations one assumes that 1 — K(f) 7 ^ 0, see |3s], 10 


More complicated is Wiener-Hopf equation 



( 21 ) 
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Assume that K(t) G L(— 00 , 00 ) and 1 — K(£) 7 ^ 0, G (— 00 , 00 ). Define 
the index: n := —— arg[l — K (O]|^°oo- The simplest result is: 

Z7T 

If and only if k — 0 equation [21\) is uniquely solvable in L p (— 00 , 00 ). 
The solution can be obtained analytically by the following scheme. De¬ 
note u + (t) i ^ _ q • Write equation (12TT) as 



K(t — s)u + (s)ds + U-{t) + / + , 


—00 < t < 00 , 


where w_ (t) = 0 for t > 0 is an unknown function. Take the Fourier transform 
of this equation to get [1— K(£)]u + = /++«_. Write 1—K(f) = A+(£)A'_(£), 
where K + and Jl_ are analytic in the upper, respectively, lower half-plane 
functions. Then 


A+K)fi + ({) = /+K)/C‘K) + fi-(O-A'-K). 

Let P + be the projection operator on the space of functions analytic in the 
upper half-plane. Then 

If u + is known, then u+(t) is known, and the Wiener-Hopf equation is solved. 
The non-trivial part of this solution is the factorization of the function 1 — 

2.2. Equations basic in random fields estimation theory 
I 11 this Subsection the results from jl 6 [ are given. 

Let u{x) = s{x) + n(x),x G M 2 ,s(x) and n{x) are random fields, s(x) is 
a “useful signal” and n(x) is “noise”. Assume that the mean values s(x) = 
n(x) = 0, the overbar denotes mean value. Assume that covariance functions 

u*(x)u(y) := R(x,y ), u*(x)s(y) := f(x,y ), 

are known. The star here denotes complex conjugate. Suppose that a linear 
estimate Lu of the observations in a domain D is 

Lu= h(x,y)u(y)dy, 

J D 
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where h(x, y ) is, in general, a distributional kernel. Then the optimal esti¬ 
mate of As, where A is a linear operator, is found by solving the optimization 
problem 

e := (Lu — As) 2 = min . (22) 

A necessary condition for the kernel (filter) h(x, y) to solve this problem is 
the integral equation: 

Rh := / R(x,y)h(z,y)dy = f(x, z), x,zED\=DUT, (23) 
J D 

where T := dD is the boundary of D. Since z enters into equation (f23l) as a 
parameter, one has to study the following integral equation 


R(x, y)h(y)dy = f(x), x e D. 


'D 


The questions of interest are: 


(24) 


1) Under what assumptions equation fl23| has a unique solution in some 
space of distributions? 

2) What is the order of singularity of this solution? 

3) Is this solution stable under small (in some sense) perturbations of /? 
Does it solve estimation problem ([22])? 

4) How does one calculate this solution analytically and numerically? 

If a distribution h = D l hi, where hi G Lj^ , then the order of singularity 
of h, ordsing/z, is |Z| if hi is not smoother than L| . ^ ^ 0116 denotes 

Q\ l \ 

h + h + • • • + lr, D l = —j-. For the distribution theory one can 

dx 1 ^ ... dx l r r 

consult, for example, [2j. 

Let us define a class of kernels R(x,y), or, which is the same, a class 1Z 
of integral equations (124|) . for which the questions 1) - 4) can be answered. 

Let L be a selfadjoint elliptic operator of order s in L 2 (M r ), A, dp( A) 
and <h(x, y, A) are, respectively, its spectrum, spectral measure and spectral 
kernel. This means that the spectral function E\ of the operator L has the 
kernel E\(x,y) = f A ^>(x,y, A)dp(A) and a function <t>(L) has the kernel 


4>{L){x,y)= [ (j>(A)$(x,y,A)dp(A), 
J A 
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where 0(A) is an arbitrary function such that f A |0(A)| 2 d(.E\/, f) < oo for 
any / G D(0(L)). 

Define class 1Z of the kernels R(x, y ) by the formula 

R{x,y) = I ^^$(x,y,\)dp(\), (25) 

where P( A) > 0 and Q( A) >0 are polynomials of degree p and, respectively, 
q, and p < q. Let a := 0.5(g — p)s. The number a > 0 is an integer because 
p and q are even integers if the polynomials P > 0 and Q > 0 are positive 
for all values of A. Denote by H a = H a (D) the Sobolev space and by H~ a 
its dual space with respect to the L 2 (D) inner product. Our basic results are 
formulated in the following theorems, obtained in (l 6 | . 

Theorem 2.1. If R(x,y) G TZ, then the operator R in equation [2f\ ) is an 
isomorphism between the spaces H~ a and H a . The solution to equation p4l ) 
of minimal order of singularity, ordsingh < a, exists, is unique, and can be 
calculated by the formula 


h(x) =Q(L)G,G(x ) = 


( 26 ) 


g(x)+v(x) in D 
u(x) in D' ■= M r \ D 

where g{x) G //°- 5 s (p+ < ?) an arbitrary fixed solution to the equation 

P(L)g = f in D, f G H a , (27) 

and the functions u(x) and v(x) are the unique solution to the problem 

Q(L)u = 0 in D', u(oo) = 0, (28) 

P(L)v = 0 in D, (29) 

dffU — & N {v + g) on dD , 0 < j < 0.5 s(p + q) — 1. (30) 

Here &nu is the derivative of u along the outer unit normal N. 

Theorem 12.11 gives answers to questions 1) - 4), except the question of 
numerical solution of equation (t24|) . This question is discussed in 


Remark 2.2. If one considers the operator R G 1Z as an operator in L 2 (D), 
R : L 2 (D) — > L 2 (D) and denotes by A j = A j(D) its eigenvalues, Ai > A 2 > 
• • ■ > 0 , then 

A j = cj~^ q ~ p ^ 3 ^ r [ 1 + o(l)] as j — > 00 , (31) 

where c = 7 ^~R> s / r , y := (27r) _r f D p(x)dx, and the rj(x) is defined in |l( 
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Example. Consider equation (EMj) with D = [—1,1], R(x,y ) = exp(— \x — 
V I), L = -i£, r = 1 , s = 1 , P{ A) = 1, Q{ A) = (A 2 + l)/2, $(x,2 /, A) = 
(27t) _ 1 exp{zA(x — y)}, and dp{ A = dA). Then 



2 / 1 )%)*/ = /, 


-1 < x < 1, 


a = 1, and Theorem 12. II yields the following formula for the solution lv. 


h(x ) = 


-/" + / , /'(!) + /(I) 


+ 


% - 1 ) + 


-/'(-!) + /(-!) 


% + 1 ), 


where d is the delta-function. 


3. Singular integral equations 

3.1. One-dimensional singular integral equations 


The results in this Section are taken from [11 
Consider the equation 


25 and 27 


Au := a(t)u(t ) H- 

*7T J L S 


M{t, s) 


u(s)ds = fit). 


(32) 


Here the functions a(f), M(t,s), and fit) are known. Assume that they 
satisfy the Holder condition. The contour L is a closed smooth curve on 
the complex plane, L is the boundary of a connected smooth domain D, 
and D' is its complement on the complex plane. Since M(t, s) is assumed 
Holder-continuous, one can transform equation (132]) to the form 


Au := a{t)u -\—— f ^ ^ + f Kit , s)u(s)ds = /(t), (33) 

vk J L s-t J L 


where the operator with kernel K(t,s ) is of Fredholm type, bit) = Ad(t,t), 
K(t, s) = _L . T^he operator A in (l33l) is of the form A = A 0 + K , 

where A 0 u = a{t)u + Su , and Su = A J L u [ b ^ b ■ The singular operator S is 

defined as the limit lim — 
e->0 m 

The operator S maps L P (L) into itself and is bounded, 1 < p < oo. Let 
us denote the space of Holder-continuous functions on L with exponent 7 by 


I" u(s)ds 

'\s—t\>e s 


t 
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Lip (L). The operator S maps the space Lip 7 (L) into itself and is bounded 
in this space. If $( 2 ) = Aq f L ai f} db and u E Lip Q (L), then 

u(s)ds 


. , . , uit) 1 

4 W =Hr + 2™ 


s — t ’ 


where +(—) denotes the limit z -E t, z E D (D 
If u E L P (L), 1 < p < 00 , then 


s/“ Whl 


s — t | 


ds = mu(t) + 


u(s)ds 
s — t 


One has S' 2 = /, that is, for 1 < p < 00 , 


(iny 


t — t 


u(s)ds 


s — T 


(34) 


dsjdr = u(t), \/u E L V (L). 


(35) 


(36) 


Let n = dim N (A) and n* := dimiV(A*). Then index of the operator A is 
defined by the formula indA n — n* if at least one of the numbers n or n* 
is finite. 

If A = B + K , where B is an isomorphism and K is compact, then A is 
Fredholm operator and indO = 0. If A is a singular integral operator (l33]h 
then its index may be not zero, 


k := indR 


27 tJ l g a(t) + b(ty 


(37) 


where arg / denotes the argument of the function /. The operator A with 
non-zero index is called a Noether operator. Formula (137|) was derived by 
F.Noether in 1921. Equation (1331) with K(x,y ) = 0 can be reduced to the 
Riemann problem. Namely, if <F(^) = A- j L j^lds, then u(t) = <F + (t) — $ _ (t) 
and $ + (t) + &~(t) = Su by formula (j5Tj) . 

Thus, equation (ITU) with K(x,y) = 0 can be written as 


a(f)($ + - $") + 6(f)($ + + $") = /, 


or as 


., , . . ait) — bit ) 

$ + (f) = ),( , + 


/ 


a(t) + b(t ) 
This problem is written usually as 


a(t) + b(t) ’ 
® + (t) = G(t)$ (t) + g(t), t E L, 


t E L. 


(38) 

(39) 
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and is a Riemann problem for finding piece-wise analytic function $( 2 ) from 
the boundary condition (J2HJ). 

Assume that a(t) ± b{t) 7 ^ 0, Vt e L, and that L is a closed smooth curve, 
the boundary of a simply connected bounded domain D. 

If one solves problem (139|) then the solution to equation (133|) (with K(x, y ) = 
0) is given by the formula u = $ + — <f> _ . Let k = indG. 

Theorem 3.1. If k > 0 then problem (j39j) is solvable for any g and its 
general solution is 


<L(z) 


X(z) r g(s)ds 
2m J L X+(s)(s - z) 


+ X(z)P K (z), 


(40) 


where 


X + (z) 

r(-) = 


= exp(T+(z)), X (2) 

j_ r 

2m J L s-z 


= z K exp(T-(z)), 


(41) 

(42) 


where P K (z) is a polynomial of degree k with arbitrary coefficients. If k = —l, 
then problem (|39li is uniquely solvable and A" - ( 00 ) = 0. If k < — 1 then 
problem (l39]h in general, does not have a solution. For its solvability, it is 
necessary and sufficient that 

[ v * k ~ l d.s = 0, k = 1, 2,..., -k - 1. (43) 

JL ^ + {S) 

If these conditions hold, then the solution to problem (f3 7 Jll is given by formula 

(HDD with p K {z ) = 0 . 

This theorem is proved, for example, in Q and 0 - It gives also an 
analytic formula for solving integral equation (l33l) with k(x,y ) = 0 . 

If one consider equation (ITIl) with k(x,y ) 7 ^ 0, then this equation can be 
transformed to a Fredholm-type equation by inverting the operator A 0 . A 
detailed theory can be found in [ 3 ] and (Toj ] . 


3.2. Multidimensional singular integral equations 
Consider the singular integral 


An : = 


r m f(x,6)u(y)dy, x,y e R n 


= \x- y I, 


(44) 


'R" 
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Assume that u(x) = O(Af) for large |x|, that k > 0, and / is a bounded 
function continuous, with respect to 6. Let f s f(x, 9)dS = 0, where S stands 
in this section for the unit sphere in R m . 

Suppose that \u(x)—u(y)\ < cir a (l+|x| 2 ) _fe//2 for r < 1, where c\ = const , 
and |w(x)| < C 2 (l + |x| 2 ) _fe / 2 , where C 2 = const. Denote this set of functions 
u by A at k. If the first inequality holds, but the second is replaced by the 
inequality |u(x)| < c 3 (l + \x\ 2 )~ k ^ 2 ln(l + |x| 2 ), then this set of u is denoted 
A' 

a,k ' 

Theorem 3.2. The operator ilTTli maps A a ^ with k < m into A' ak , and A a ^ 
with k > m into A a)7n . 

Let D C R m be a domain, possibly D = R m , the function /(x, 6) be 
continuously differentiable with respect to both variables, and u satisfies 
the Dini condition, that is sup|| a ,_ y || <( 5 ;I , jygi: , \u(x) — u(y)\ = w(u,t), where 

fo —dt < oo, 5 = const > 0. 

Then 


d f f(x,9) 

—u[y)dy = 


dxj J D rm 


Assume that 


[ d u{v) £ 


3 L 


f(x,0) 

r m—l 


dy-u{x) I f(x,0)cos(r Xj )ds. 

(45) 


| f (x, 9)\ p, dS <C 0 = const , - + - = 1, 

p p' 


( 46 ) 


and let 

A 0 u = f — V jP~u{y)dy. ( 47 ) 

JR” r 

Theorem 3.3. If (1451) holds then 


IIAjIIlp < csu pII/(^,^)IIlp'( S )- 

x 

The operator (1471) is a particular case of a pseudodifferential operator 
defined by the formula 

Au= f f e l{x ~ y)< cr A (x,f)u(y)dydf, (48) 

J R m J R m 

where the function cr^x, £) is called the symbol of A. For the operator 
Au = a(x)u(x) + A 0 u, where A 0 is defined in formula (147)) . the symbol of A 
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is defined as a(x,£) = a(x) + K(x, £), where K(x,x — y ) = and K( x,£) 
is the Fourier transform of K(x,z) with respect to z. One has 


cta+b{x, 0 = a A (x, f) + a B (x, £), a AB (x , f) = a A (x, £)(t b {x, £)■ 


It is possible to estimate the norm of the operator (H7|) in terms of its 
symbol, 

<7Ao(x,£) = f ^ e~ iy ' z dy, y° = y- £ = 

J r- \y\ m \y I N 

If T is a compact operator in L 2 (R m ) then its symbol is equal to zero. Let 
Au = a{x)u + A 0 u, and a A (x, £) = a(x) -\-cr Ao (x, ^). If a A (x,£,) is sufficiently 
smooth and does not vanish, then the singular integral equation Au + Tu = f 
is of Fredholm type and its index is zero. In general, index of a system 
of singular integral equations (and of more general systems) is calculated 
in the work [l|. Additional material about singular integral equations and 


pseudo differential operators one finds in [11] and 


4. Nonlinear integral equations 

The results from this Subsection are taken from 
Let us call the operators 


3 , 0 . 0 - 


Uu — / K(x,t)u(t)dt, Hu= / K(x,t)f(t,u(t))dt , 

J d J D 

and Fu := f(t,u(t )), respectively, Urysohn, Hammerstein, and Nemytskii 
operators. 

The operator F acts from L p := L P (D) into L q := L q (D) if and only if 

\f(t,u)\ < a(t) +b\u\ p/q , p < oo, a(t) £ L q , b = const. 

If p — oo, then F : L°° —> L q if and only if \f(t,u)\ < cih(t), a,h(t ) G L q , 
and \u(t)\ < h, 0 < h < oo. The function f(t,u), t £ D, u G (—oo,oo), 
is assumed to satisfy the Caratheodory conditions, that is, this function is 
continuous with respect to u for almost all t £ D and is measurable with 
respect to t £ D for all u £ (— 00 , 00 ). If F : L p ^ L q and f(t,u ) satisfies 
the Caratheodory conditions, then F is compact provided that one of the 
following conditions holds: 1) q < 00 , 2) q = 00 , p < 00 , f(t,u ) = a(t), 3) 
q = p = 00 , | f(t, u ) — f(t, u)| < 4>h(u — v ) with respect to z for any h > 0. 
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If / satisfies the Caratheodory conditions, and K(x,t ) is a measurable 
function on D x D, and F : L p —>■ L r , K : L r —> L q , then the operator 
F[ = KF : L p — y L q , F[ is continuous if r < oo, and F[ is compact if r < oo 
and K is compact. These results one can find, for example, in 27 . 

A nonlinear operator A : X —> Y is Frechet differentiable if 


A(u + h) — A(u ) = Bh + w(u, h ), 

where B : X —> Y is a linear bounded operator and lining | x _> 0 
Consider the equation 


\\w(u,h) ||y _ g 


u(x) — n / K(x,t,u(t))dt + f(x) := Au, (49) 

J D 

where D C R m is a bounded closed set, meas-D > 0, the functions K(x,t,u) 
and f(x) are given, the function u is unknown, and fi is a number. 

Theorem 4.1. Let us assume that A is a contraction on a set M, that is, 
||.A(w) — A(u) < q\\u — n||, where u, v G M, M is a subset of a Banach space 
X, 0 < q < 1 is a number, and A : M —>■ M. Then equation (1491) has a 
unique solution u in M, u = lim^ooU^ , where u n+ \ = A(u n ),uo G M, and 
\\u n - u\\ < j^\\ui - u 0 \\,n = 1 , 2 ,.... 

This result is known as the contraction mapping principle. 

Theorem 4.2. Assume that A is a compact operator, M <Z X is a bounded 
closed convex set, and A : M M. Then equation (149]) has a solution. 

This result is called the Schauder principle. The solution is a fixed point 
of the mapping A, that is, u = A(u). Uniqueness of the solution in Theorem 
14.21 is not claimed: in general, there can be more than one solution. 

Remark 4.3. A version of Theorem 14.21 can be formulated as follows: 

Theorem 4.4. Assume that A is a continuous operator which maps a convex 
closed set M into its compact subset. Then A has a fixed point in M. 

Let us formulate the Leray-Schauder principle: 

Theorem 4.5. Let A be a compact operator, A : X —» X, X be a Banach 
space. Let all the solutions u(X) of the equation u = A{u,X), 0 < A < 1, 
satisfy the estimate sup 0<A<1 ||w(A)|| < a < oo. Let ||A(u,0)|| < b for b > a 
and ||u|| = b. Then the equation u = A(u, 1) has a fixed point in the ball 
| \u \| < a. 
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Let Del be a bounded convex domain, A n u C D if u G D, n = 1, 2,.. 

D is the closure of D , A n rt 7 ^ u for n > no and u G dD, where dD is the 
boundary of D. Finally, let us assume that A is compact. 


Theorem 4.6. Under the above assumptions A has a fixed point in D. 


This theorem is proved in [7] 

Consider the equation u = H{u),H : L 2 (D ) —> L 2 (D), and assume that 
K(x,t) is positive-definite, continuous kernel, f(t,u ) is continuous with re¬ 


spect tot E D and u G (— 00 , 00 ) function, such that f Q u f(t, s)ds < 0.5 au 2 +b, 
where a < A, and A is the maximal eigenvalue of the linear operator K in 
L 2 (D). Then there exists a fixed point of the operator H. This result can 


be found in 27 


5. Applications 


13 , 221 and 24 


5.1. General remarks. 

The results in this Subsection are taken from 

There are many applications of integral equations. In this section, there 
is no space to describe in detail applications to solving boundary value prob¬ 
lems by integral equations involving potentials of single and double layers, 
applications of integral equations in the elasticity theory ([lo|, [9j]), applica¬ 
tions to acoustics electrodynamics, etc. 

We will restrict this Subsection to some questions not considered by other 
authors. The first question deals with the possibility to express potentials of 
the single layer by potentials of the double layer and vice versa. The results 
are taken from [13] . The second application deals with the wave scattering 
by small bodies of an arbitrary shape. These results are taken from 
24 . 


221 and 


5.2. Potentials of a single and double layers. 

r gik\x—t\ 

V(x)= g(x,t)a(t)dt, g(x,t) = (50) 

Js An\x-t\ 

where S' is a boundary of a smooth bounded domain D C with the 
boundary S, a(t) G Lip 7 (S), 0 < 7 < 1. 

It is known that V(x) G C'(R 3 ), V(oo) = 0, 


V£(a) = 


Aa ± a(s) 


Aa := 2 


dg(s,t) 

dN s 


a(t)dt, 


(51) 
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where N = N s is the unit normal to S at the point s pointing out of D, +(—) 
denotes the limiting value of the normal derivative of V when x —> s G S 
and x G D (£>'), where D' := R 3 \ D. 

Potential of double layer is defined as follows: 



One has the following properties of W: 



(52) 

(53) 


Theorem 5.1. For any V(W) there exists a W(V) such that W = V in D. 
The V{W) is uniquely defined. 

Theorem 5.2. A necessary and sufficient condition for V(cr) = W([f) in D' 
is 

/ Vhjdt = 0, 1 < j < r', (/ + A)hj = 0, (54) 

Js 

where the set {hj} forms a basis of N(I + A), and A is defined in (1511) . 

A necessary and sufficient condition for W (/i) = V (cr) in D' is: 

/ Wujdt = 0, 1 < j < r, Aaj — <7j = 0, (55) 

Js 

where the set {af} forms a basis of N (/ — A). 

Theorem 15.11 and 15.21 are proved in [l^l . 

5.3. Wave scattering by small bodies of an arbitrary shape. 

The results in this Subsection are taken from j22| and [24|. Consider the 
wave scattering problem: 


(V 2 + k 2 )u = 0 in D' := R 3 \ D, = (u on S = dD, 
u = u 0 + v, u 0 = e lka ' x , a G S 2 , 


(56) 

(57) 

(58) 



r = he —> oo. 
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Here k = const > 0, S 2 is the unit sphere in R 3 , ( = const is a given 
parameter, the boundary impedance, D is a small body, a particle, S is its 
boundary, which we assume Holder-continuous, and N is the unit normal to 
S pointing out of D. If Imf < 0 then problem (l56]) - (l58j) has exactly one 
solution, v is the scattered held, 

e ikr /f\ x 

v = A(/3, a , k) -h o ( - ) , HI = r —y oo, f3 — —, (59) 

r \r J r 


A(j3, a , k) is called the scattering amplitude. 

Our basic assumption is the smallness of the body D: this body is small 
if ka -C 1, where a = 0.5 diamD. 

In applications a is called the characteristic size of D. 

Let us look for the solution to problem (|56l) - (|58l) of the form 


r gik\x—1\ 

u(x) = U 0 (x) + / g(x,t)a(t)dt, g(x,t) = — --, 

Js An\x-t\ 


and write 


(60) 


g(x,t)a(t)dt = g(x,x,)Q + / \g(x,t) - g(x,X!)]cr(t)dt, (61) 
Js Js 

where X\ G D is an arbitrary fixed point inside D, and 


Q'= cr(t)dt. 

Js 


If ka <C 1 and \x — X\\ := d a, then 


(62) 


\g(x,t) - g(x,Xi)}cr(t)dt 


« \g(x,Xi)\Q. 


(63) 


Therefore the scattering problem (l56l) - (l58l) has an approximate solution of 
the form: 


u(x) — uq(x) + g(x,xi)Q, ka 1, \x — xf\^>a. (64) 

The solution is reduced to finding just one number Q in contrast with the 
usual methods, based on the boundary integral equation for the unknown 
function aft). 
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To find the main term of the asymptotic of Q as a —$■ 0, one uses the 
exact boundary integral equation: 

uqn{s ) - C«oO) + Aa ~ C f 9 (s, t)a(t)dt = 0, s G S, (65) 

J s 

where formula (15T1) was used. 

We do not want to solve equation (1551) which is only numerically possible, 
but want to derive an asymptotically exact analytic formula for Q as a —> 
0. Integrate both sides of formula (155]) over S. The first term is equal to 
f s u 0N ds = f D V 2 u 0 dx ~ 0(a 3 ). The sign ~ stands for the equality up to 
the terms of higher order of smallness as a —» 0. The second term is equal 
to —( f s u 0 ds ~ —Cm 0 (xi)| j S'|, where |,S| = 0(a 2 ) is the surface area of S. 
The third term is equal to \ f s Aads — \Q. When a —> 0 one checks that 
Aa ~ A 0 a, where A 0 = ^4|fc=o, and f s A 0 adt — — f s adt. Therefore the third 
term is equal to — Q up to the term of higher order of smallness. The fourth 
term is equal to — ( J s dta(t ) J s g(s,t)ds = o(Q), as a —* 0. Thus, 

Q - -C|*S1no(ar), a ->• 0, (66) 

and Mo(^i) — 1, 

Af o \ C|5'K(a:i) 

a, K) = ---• (67) 

47T 

For the scattering by a single small body one can choose X\ G D to be the 
origin, and take Mo(^i) = 1. But in the many-body scattering problem the 
role of Wo is played by the effective held which depends on x. Formulas 
(154]) -f??f solve the scattering problem (156]) - (158|) for one small body D of an 
arbitrary shape if the impedance boundary condition (156]) is imposed. The 
scattering in this case is isotropic and A(/3, a, n) = 0(a 2 (). 

If the boundary condition is the Dirichlet one, u\s = 0, then one derives 
that 

Q ~ — Cuo, a —> 0; wo — 1, (68) 

and 

C 

A(/3,a,K) = -, a —» 0, u 0 ~ 1, C = 0(a ). (69) 

47T 

Here C is the electrical capacitance of the perfect conductor with the shape 
D. 

In this case the scattering is isotropic and |H| = O(a), a —$■ 0. Thus, the 
scattered held is much larger than for the impedance boundary condition. 
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For the Neumann boundary condition % = 0 on S the scattering is 
anisotropic and A(f3, a, k) = 0(a 3 ), which is much smaller than for impedance 
boundary condition. 

One has for the Neumann boundary condition the following formula: 

A(P, a, k) = ^ (ik/3 pq ^-/3 p + V 2 u 0 ^ , (70) 

where f3 pq is some tensor, \D\ is the volume of \D\, (3 p = lim^i-nx, —, x p := 
x ■ e p is the p— th coordinate of a vector x G R 3 , u 0 and its derivatives are 
calculated at an arbitrary point inside D. This point can be chosen as the 
origin of the coordinate system. Since D is small, the choice of this point 
does not influence the results. The tensor /3 pq is defined for a body with 
volume V and boundary S as follows: 

1 f 

Ppq = — J t p CT q dt, 

where a q solves the equation a q = Aa q — 2 N q , N q N ■ e q , and {e q } 3 =1 is a 
Cartesian basis of R 3 . 

Finally, consider the many-body scattering problem. Its statement can 
be written also as (|56 ]) - (|58j) but now D = U ^ =1 D m is the union of many 
small bodies, C — Cm on S rn — dD m . 

One looks for the solution of the form 


M 


u(x) = u 0 (x) + g(x,t)a m (t)dt. 


m =1 ' 


(71) 


Let us dehne the effective held in the medium by the formula 


u e (x) = u Q {x) + ^g(x,x m )Qm, \x-x m \~a. (72) 

m^j 

Here x m G D m are arbitrary fixed points, and \x — x m \ ~ a means that 
\x — x m \ is of the size of a. 

Assume that these points are distributed by the formula: 

A/~(A) = — f N(x)dx[l + o(l)], a — > 0. (73) 

o K Ja 
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Here AT (A) = £ ImGA 1 is the number of points in an arbitrary open set A, 
N(x) > 0 is a given function, 0 < k < 1 is a parameter, and an experimen¬ 
talist can choose N(x) and k as he wishes. Let us assume that ( m = , 

where h(x) is an arbitrary continuous in D function such that Imh(x ) < 0. 
This function can be chosen by an experimenter as he wishes. 

Under these assumptions in [24j the following results are proved. 

Assume that \S m \ — ca 2 , where c > 0 is a constant depending on the shape 
of D m . Then Q m ~ -ch(x m )u e (x m )a 2 ~ K . Denote u e (x m ) := u m ,h(x m ) := 
h m . Then one can find the unknown numbers u m from the following linear 
algebraic system (LAS): 


M 

Uj U(yj C ^ ^ gjm,hmU m a 
m+j 


1 <j<M, 


9 jm ■ 


e ik\xj-Xm\ 

47 t\xj - x m 


(74) 


LAS (l74l) can be reduced to a LAS of much smaller order. Namely, denote by 
D a finite domain in which all the small bodies are located. Partition D into 
a union of P cubes A p with the side b = b(a), the cubes are non-intersecting 
in the sense that they do not have common interior points, they can have 
only pieces of common boundary. Denote by d = d(a) the smallest distance 
between neighboring bodies D m . Assume that 


a < d C 6, ka 1. (75) 

Then system (1731) can be reduced to the following LAS: 

p 

Uq Uq q ^ fjqpkp llpAp | Ap |, 1 A Q A P) (76) 

where N p := N(x p ), x p G A p , is an arbitrary point, N(x) is the function from 
formula (173]) . and |A p | is the volume of the cube A p . 

The effective field has a limit as a —> 0, and this limit solves the equation 


u(x ) = u 0 (x) -c g(x,y)h(y)N(y)u(y)dy, (77) 

J n 

where the constant c is the constant in the definition |,S m | = ca 2 . One may 
consider by the same method the case of small bodies of various sizes. In this 
case c = c m , but we do not go into details. If the small bodies are spheres of 
radius a then c = An. 
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Applying the operator V 2 + k 2 to the equation (1771) one gets 

(V 2 + k 2 n 2 (x))u := (V 2 + k 2 — cN(y)h(y))u = 0 in R 3 . (78) 

Therefore, embedding of many small impedance particles, distributed accord¬ 
ing formula (1731) leads to a medium with a zero refraction coefficient. 

n 2 (x) = 1 — k~ 2 cN (x)hfx). (79) 

Since N(x) and h(x) can be chosen as one wishes, with the only restrictions 
N(x) > 0, Imh(x) < 0, one can create a medium with a desired refraction 
coefficient by formula (179|) if one chooses suitable N(x) and h(x). 

A similar theory is developed for electromagnetic wave scattering by small 
impedance bodies in j23|. 
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